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O 1 ABSTRACT 

\Jj Aims. Brown dwarfs are covered by dust cloud layers which cause inhomogeneous surface features and move below the observable r = 1 level 
l ^ l during the object's evolution. In any case, the cloud layers have a strong influence on the structure and spectral appearance of brown dwarfs 
and extra-solar planets, e. g. by providing large local opacities and by removing condensable elements from the atmosphere causing a sub-solar 
< metalicity in the atmosphere. We aim at understanding the formation of cloud layers in quasi-static substellar atmospheres which consist of 
^ ' dirty grains composed of numerous small islands of different solid condensates. 

t Methods. The time-dependent description presented in (Helling & Woitke 2006) is a kinetic model describing nucleation, growth and evapora- 
• tion. It is extended to treat gravitational settling and is applied to the static-stationary case of substellar model atmospheres. From the solution 
for the dust moments, we determine the grain size distribution function approximately which, together with the calculated material volume 
X J # , fractions, provides the basis for applying effective medium theory and Mie theory to calculate the opacities of the composite dust grains. 

1 Results. The cloud particles in brown dwarfs and hot giant-gas planets are found to be small in the high atmospheric layers (a ss 0.01 yum), and 
composed of a rich mixture of all considered condensates, in particular the abundant MgSiO^ [s], Mg2Si04[s] and Si02[s]. As the particles settle 
OO , downward, they increase in size and reach several 100 yum in the deepest layers. The more volatile parts of the grains evaporate and the particles 
■ stepwise purify to form composite particles of high-temperature condensates in the deeper layers, mainly Fe[s] and A^Osfs]. The gas phase 
^ \ abundances of the elements involved in the dust formation process vary by orders of magnitudes throughout the atmosphere. The grain size 
• <"H , distribution is found to be relatively broad in the upper atmospheric layers but strongly peaked in the deeper layers. This reflects the cessation 
/\l ' of the nucleation process at intermediate heights. The spectral appearance of the cloud layers in the mid IR (7 - 20yum) is close to a grey 
^ body with only weak broad features on a few percent level, mainly caused by MgSiOsfs], and Mg2SiO,i[s]. These features are, nevertheless, a 
- - fingerprint of the dust in the higher atmospheric layers that can be probed by observations. 

Conclusions. Our models predict that the gas phase depletion is much weaker as compared to phase-equilibrium calculations in the high atmo- 
spheric layers. Because of the low densities, the dust formation process is incomplete there, which results in considerable amounts of left-over 
elements that might produce stronger and broader neutral metallic lines. 

Key words. Stars: atmospheres, Stars: low mass, brown dwarfs, methods: numerical, astrochemistry 



1. Introduction 

Dust in form of small solid particles (grains) becomes an in- 
creasingly important component in understanding the nature 
of substellar objects with decreasing T^, i.e. brown dwarfs 
and giant-gas planets. Observations start to provide direct ev- 
idence for dust clouds covering brown dwarfs (Cushing et al. 
2006) and extrasolar giant-gas planets (Richardson et al. 2007, 
Swain et al. 2007). The search for biosignatures in extraterres- 
trial planets becomes more complicated if such cloud layers 
cover the atmosphere that efficiently absorb in the wavelength 
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region where e.g. the Earth vegetation's red edge spectroscopic 
features (600 - 1 100 nm, Saeger et al. 2007) or the extraterres- 
trial equivalents are situated. In fact, the development of live 
seems impossible below optically thick cloud layers, because 
the star light needs to reach the surface to create the neces- 
sary departures from thermodynamical equilibrium that allow 
for structure formation. Furthermore, abundances of molecules 
like O2 and O3 as the carriers of spectral biosignatures will 
be strongly affected by the presence of cloud layers, because 
of chemical surface reactions and the element depletion due to 
dust formation in the atmosphere. Thus, the understanding of 
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the details of cloud formation physics and chemistry is a major 
issue in modelling substellar atmospheres. 

This paper presents a kinetic approach for modelling quasi- 
static atmospheres with stationary dust cloud layers includ- 
ing seed particle formation (nucleation), grain growth, grav- 
itational settling, grain evaporation, and element conserva- 
tion (Sects. |2] and [3). The model is a further development of 
the time-dependent description presented in Helling & Woitke 
(2006; Paper V) and is particularly suited to treat the forma- 
tion of "dirty" dust grains (i. e. particles composed of numer- 
ous small islands of different solid condensates) in the frame of 
classical stellar atmospheres with consistent radiative transfer 
and convection (Dehn 2007, Helling et al. 2008). The results 
of the models (see Sect. [4} describe the vertical cloud structure 
and the amount of dust formed in the atmosphere of substellar 
objects as well as the amount of condensable elements left in 
the gas phase. The models provide further details like the ma- 
terial composition of the cloud particles, the mean grain sizes 
and the size distributions as function of atmospheric height. 
Section |5] demonstrates what spectral features such cloud lay- 
ers made of dirty solid particles exhibit, and from which tem- 
perature and pressure level these features originate. 

2. Nucleation, growth and evaporation of 
precipitating dirty grains 

In Paper V, a kinetic description for the nucleation, growth 
and evaporation of dirty dust particles has been devel- 
oped by extending the time-dependent moment method of 
Gail & Sedlmayr (1988). The basic idea is that a "dirty" solid 
mantle will grow on top of the seed particles, because these 
seeds can only form at relatively low temperates (T 5 1400 K) 
where the oxygen-rich gas is strongly supersaturated with re- 
spect to several solid materials. The dirty mantle is assumed to 
be composed of numerous small islands of different pure con- 
densates. The formation of islands is motivated by experiments 
in solid state physics (Ledieu et al. 2005, also Zinke-Allmang 
1999) and by observations of coated terrestrial dust particles 
(Levin et al. 1996, Korhonen et al. 2003). Note that we con- 
sider dust formation by gas-solid reactions only and omit solid- 
solid reactions and lattic rearrangements inside the grains in 
our model. 

In the following, we extent this description to include the 
effects of gravitational settling (drift, rain-out, precipitation) 
which is important to understand the long-term quasi-static 
structures of brown dwarf and gas-giant atmospheres. The chal- 
lenge here is to account properly for the element conservation 
when dust particles consume certain elements from the gas 
phase at the sites of their formation, transport them via drift 
motions through the gas and finally release the elements by 
evaporation at other places. Furthermore, we want to abandon 
the assumption that the number of elements equals the number 
of condensates made in Paper V, because there are typically 
much more condensates than elements. 

We consider the moments Lj(x, t) [cm J g _1 ] (j = 0, 1, 2, ...) 
of the dust volume distribution function f(V, x, t) [cm -6 ] (for 
more details, see Paper V, Sect. 2.1), where x and t are space 
and time. V [cm 3 ] is the volume of an individual dust particle. 



The total dust volume per cm 3 stellar matter, V tot , is given by 
the 3 rd dust moment as 

P L 3 = f(V)VdV = V tot [cm 3 cm -3 ], (1) 

where p [gcrrT 3 ] is the mass density and V( is the lower inte- 
gration boundary. In an analogous way, we define the volume 
V s of a certain solid species s by 

pL" 3 = I f(V)V s dV = V s [cm 3 cnT 3 ], (2) 

JV[ 

where V s [cm 3 ] is the sum of island volumes of material s 
in one individual dust particle. For simplicity, we assume that 
V s /V = Vs/Vt t is constant for all dust particles at a certain 
position in the atmosphere, i. e. we assume a unique volume 
composition of all grains at (x, f), such that V tot = 2 V and 

By means of this assumption, it is possible to express the in- 
tegrals that occur after integrating the master equation (Eq. (1) 
in Paper V) over size in terms of other moments. The results 
are the dust moment conservation equations. The change of 
the partial volume of solid s can then be expressed analog to 
Eq. (23) in Paper V 

J t (pL s 3 ) + v(j^7(V)V s (v gas +v^V))rfv) 

= V//* +xL,pLi (3) 

where we have assumed that the Knudsen numbers are large 
(Kn » 1, compare Paper II) and that the drift velocities Vd r ( V) 
of the dust particles can be approximated by the equilibrium 
drift velocities v di .(V) (final fall speeds). 

The source terms on the r.h.s. of Eq. (f3]) describe the ef- 
fects of nucleation, growth and evaporation of condensate s. 

= J(V ( ) = f(Ve)^\ v=v [s-'cirT 3 ] is the stationary nucle- 
ation rate (see Sect. [33). Ve s [cm 3 ] is the volume occupied 
by condensate s in the seed particles when they enter the in- 
tegration domain in size space. The net growth velocity of con- 
densate s Xnet [cms -1 ] (negative for evaporation) is given by 
(Eq.(24) in Paper V) 

X nel =V36n)^ (4) 

r =l V r \ ^"surf/ 

Here, r is an index for the chemical surface reactions (see 
Table [TJ, A V r s is the volume increment of solid s by reaction 
r (Yj AV, s = AV r ), «r ey is the particle density of the key reac- 
tant, v^ 1 is its thermal relative velocity and a r is the sticking 
coefficient of reaction r. S r is the reaction supersaturation ratio 
and Z? s s urf = V to t/V s is a ^-factor which describes the probabil- 
ity to find a surface of kind s on the total surface. Putting Z? s s urf 
independent of V, we assume that all grains at a certain point 
in the atmosphere have the same surface and volume composi- 
tion, i.e. the grain material is a homogeneous mix of islands of 
different kinds (for more details see Paper V). 

The divergence of the drift term in Eq. (0 is treated in the 
following way. Assuming large Knudsen numbers and sub- 
sonic drift velocities, the equilibrium drift velocity is given 
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by v dr (V) = -e z ay/jrgp A /(2pcr) (Eq. 63 in Paper II), where 
a is the particle radius, e z the vertical unit vector (pointing 
upwards) and g the gravitational acceleration (downwards). 
Pi = ZPs Vs/Vtot ' s t ' le dirty dust material density and p s the 
material density of a pure condensate s. cj — sj2kT/p is a mean 
thermal velocity with T being the temperature, k the Boltzmann 
constant and p the mean molecular weight of the gas particles. 
Inserting this formula into the drift term in Eq. © yields 

f(V) V s v dl .(V) dV = - e z fa — tt^U (5) 

JV[ C T ^tot 

with the abbreviation £iK n = {j^f^ ~yS (note the difference to 
Eq.(3) in Paper III). Defining^ = L A V s /V tot , Eq.© in this 
paper can be rewritten as 

j t (pLl)+V(pI% v gas ) = VM+x^pLi + fiKn£ f^) (6) 

Equation © describes the evolution of the partial dust vol- 
ume of solid s in space and time due to advection, nucleation, 
growth, evaporation and drift. 

The third moment equation for the total dust volume pLj, 
(Eq. (1) in Paper III) can be retrieved by summing up the con- 
tributions from all condensates s as given by Eq. ©, because 

s s s s 

2.1. Quasi-static, stationary case 

In the case of a plane-parallel quasi-static stellar atmosphere 
Vg as = and the dust component is stationary j^(pLj) = 0, i. e. 
the l.h.s.of the moment equations vanish. Introducing a convec- 
tive mixing on time scale we have derived the following 
equations in Paper III for this case (see Eq. 7 in Paper III) 

dz \C T / £lKn \ T mix 3 j 

These are the moment equations with respect to the total dust 
volume for dirty grains and we use them for j — 0, 1 , 2. As out- 
lined in (Eq. (9) in Paper III), t^x is the timescale for mix- 
ing due to convective motions and overshoot which decreases 
rapidly above the convective layers with increasing height in 
the atmosphere. However, instead of using Eq. ([8]) with j = 3 
for one pure condensate, we have to use a set of equations for 
dirty grains, i.e. the third dust moment equations for all volume 

1 Our mixing approach is very much simplified, assuming that 
the gas/dust mix at height z is exchanged by dust-free gas from the 
deep interior of the object with element abundances e° on a mixing 
timescale r mix (z) which is adapted to the results of 3D hydrodynamical 
models (Ludwig et al. 2003, 2006). Other works consider a diffusive 
mixing here (Ackerman & Marley (2001) in Cushing et al. (2007); 
Rossow (1978) in Warren et al. (2007)). Interestingly, our model ap- 
proach and the approach by Ackerman & Marley (2001) yield very 
comparably results f.i. regarding the location of the cloud layer and the 
maximum dust-to-gas-ratio as can be seen from a comparative study 
of dust cloud models (Helling et al. 2008b). 



contributions, that is one equation for each condensate s taken 
into account. From (Eq© we find 

dZ \C T / §kn \ T mix / 

Equations ® for j e {0, 1,2} and Eqs. © for s e {1,2, ...,S} 
(S is the number of solid condensates taken into account) form 
a system of (S+3) ordinary differential equations for the un- 
knowns \L\,L,2, L3, L s 4 }. 

The element conservation equations are not affected by the 
drift motion of the dust grains. Therefore, Eq. (29) in Paper V 
remains valid, from which we derived in the static stationary 
case (compare Eq. (8) in Paper III) 

"(H)(gf ~ g) . 

T mix R key r el / 1 \ 

where i enumerates the elements. Nc is the number of 
monomers in the seed particles when they enter the size in- 
tegration domain and v, is the stoichiometric coefficient of the 
seeds (TiC>2 seeds: 1 for i = Ti and 2 for i = O). v, i5 is the stoi- 
chiometric coefficient of element i in solid material s. 

The element conservation equations (Eq. ITOb provide alge- 
braic auxiliary conditions for the ODE system (Eqs. [8]& [9jl in 
the static stationary case, i. e. one first has to solve the system of 
non-linear algebraic Eqs. ( TTOb for e, at given {La,L s 4 } (the dust 
volume composition Z? s s urf is known from L s 4 ) before the r.h.s. 
of the ODE-equations can be calculated. Since J+, nf 1 and in 
particular S r , however, depend strongly on e,, this requires a 
complicated iterative procedure which creates the most prob- 
lems in practise. 

2.2. Grain size distribution function 

The dust opacity calculations with effective medium theory and 
Mie theory (see Sect© require the dust particle size distri- 
bution function f(a) [cm -4 ] at every depth in the atmosphere, 
where a [cm] is the particle radius. This function is not a direct 
result of the dust moment method applied in this paper. Only 
the total dust particle number density n^ - pL^ and the mean 
particle size (a) = -^3/(4-7:) Li/Lq are direct results that have 
been used in Helling et al. (2006). In this paper, we reconstruct 
f(a) from the calculated dust moments Lj (J = 1 ... 4) in an ap- 
proximate way. We want to avoid the zeroth moment, because 
it is only determined by a closure condition. The idea is to in- 
troduce a suitable functional formula for f(a) with a set of four 
free coefficients, and then determine these coefficients from the 
known dust moments. For further details, see AppendixlAl Two 
possible functions for f(a) are discussed in Appendix [AT] and 
IA21 

2.3. Closure condition 

Since Lq appears only on the r.h.s. of Eq. (JHJ for j-Q, we need 
a closure condition in the form Lq = Lq{L\ , L2, L3, L4) to solve 
our ODE-system. In this paper, we use the results of the size 
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distribution reconstruction technique explained in Appendix lAl 
in application to the double delta-peaked size distribution func- 
tion and write the zeroth dust moment as 

pL =N l+ N 2 . (11) 

where the two dust particle densities N\ and /V 2 [cm -3 ] are in- 
troduced in Appendix lA.il 

3. Material equations and input data 

To solve our model equations, we need an atmospheric 
(T, p, v conv )-structure and additional material equations to cal- 
culate the number densities of the key reactants n^ v , the nu- 
cleation rate 7*, the reaction supersaturation ratios S r , and the 
sticking probabilities a r . 

3.1. Atmospheric (T, p)-structure 

The approach of this paper is to investigate the behaviour of the 
dust component in quasi-static substellar atmospheres and to 
propose how our approach to treat the micro-physical processes 
of the formation of composite dust grains (see Sec. [2j can be 
included in the frame of stellar atmosphere codes. 

For this purpose, it is sufficient to use a given (T, p, v conv ) 
structure representing a brown dwarf or giant-gas atmosphere. 
The atmospheric structure used for our models are cond AMES 
atmosphere structures^. The feedback of the dust formation and 
the presence of the dust on the atmospheric structure is thereby 
neglected in this paper, although the model code to calculate 
the (T,p, Vconv) structures did incorporate some dust modelling 
(Allard et al. 2001). Such a decoupled approach may not be en- 
tirely satisfying, but the understanding of the physics of cloud 
layers requires some basic studies, before the feedback mech- 
anisms can be attacked in the frame of highly nonlinear stellar 
atmosphere codes. For consistent models of our dust treatment 
in the frame of the PHOENIX stellar atmosphere code includ- 
ing radiative transfer and convection see ( Dehn 2007, Helling 
et al. 2008a). 

3.2. Gas-phase chemistry 

We calculate the particle densities of all gaseous species, in- 
cluding n^ ey as described in Paper III according to pressure, 
temperature and the calculated, depth-dependent element abun- 
dances e, in chemical equilibrium. For the well-mixed, deep el- 
ement abundances e? we use solar abundances according to the 
cond AMES input model. For those elements that are not in- 
cluded in the calculations (abundances are calculated for Mg, 
Si, Ti, O, Fe, Al, Ca, S), we put e, = ef. 

3.3. Nucleation 

The nucleation rate is calculated for (TiO^-clusters ac- 
cording to Eq. (34) in Paper V, applying the modified classical 
nucleation theory of Gail et al. (1984). We use the value of the 

2 ftp.ens-lyon.fr/pub/users/CRAL/fallard/ 



surface tension cr fitted to small cluster data by Jeong (2000) as 
outlined in Paper III. 

3.4. Growth/Evaporation: choice of solid material and 
condensation experiments 

Several dozens of solid species have been treated in phase 
equilibrium in the literature (Sharp & Huebner 1990, Fegley 
& Lodders 1994). We have taken into account 12 solids 
(Ti0 2 [s], Al 2 3 [s], CaTi0 3 [s], Fe[s], FeO[s], FeS[s], Fe 2 3 [s], 
SiO[s], Si0 2 [s], MgO[s], MgSi0 3 [s], Mg 2 Si0 4 [s]) in phase- 
non-equilibrium to calculate the formation and composition of 
the dirty grains. Our selection is guided by the most stable con- 
densates which yet have simple stoichiometric ratios that en- 
sure that these solids can be easily built up from the gas phase. 
The selection covers the main element sinks during dust for- 
mation. Our choice is furthermore inspired by the following 
experiments. 

SiQ 2 [s], MgO[s], FeO[s], Fe 2 Q 3 [s], MgSiQ 3 [s], Mg 2 SiQ 4 [s]: 
Vapour phase condensation experiments offer strong arguments 
against equilibrium assemblages even under controlled condi- 
tions in the terrestrial laboratory. Rietmeijer et al. (1999) have 
shown that from a Fe-Mg-SiO-H 2 vapour only Mg x Fe y O z - 
condensates with simple stoichiometric ratios form during the 
condensation process in the laboratory, which is in accordance 
with some simple phase equilibrium calculations (e.g. Sharp 
& Huebner 1990), but is in contrast to well-established text 
books like (Lewis 1997). Rietmeijer et al. (1999) report on 
the formation of the end-member oxides FeO[s], Fe 2 3 [s], 
MgO[s] and Si0 2 [s], and the fosterite olivine Mg 2 Si0 4 [s]. The 
appearance of Si0 2 [s] shows that kinetic factors can favour the 
formation of meta- stable states which is in contrast to phase 
equilibrium calculations. 

Ferguson & Nuth (2006) revised the vapour pressure of 
SiO[s] and discuss implications for possible (SiO)N nucleation. 
However, the experiments of Rietmeijer et al. (1999) which 
utilise an Fe-Mg-SiO-H 2 vapour do not show any formation 
of SiO[s]. John (2002) pointed out that SiO[s] formation may 
proceed via Si0 2 [s] + Si[s] — > 2 SiO[s]. Since we are not yet 
able to treat solid-solid reactions, we must omit these processes 
in our model. 

Fe[s], FeS[s]: We add Fe[s] as high temperature condensate 
and FeS[s] as possible sulphur binding solid in accordance with 
the experimental findings by Kern et al. (1993) and Lauretta et 
al. (1996). 

Ti0 2 [s], Al 2 3 [s], CaTi0 3 [s]: (T10 2 )n acts as nucleation 
species in our model (seed particle formation). For consistency, 
Ti0 2 [s] must therefore also be included as high-temperature 
solid compound to account correctly for the depletion of Ti 
from the gas phase. CaTi0 3 [s] is a very stable condensate 
and is included to study the consumption of Ca from the gas 
phase. According to our knowledge, no laboratory measure- 
ments on Ca-Ti oxides are available so far. The experiments by 
Kern et al. (1993) suggest the formation of Al 2 3 [s] as high- 
temperature condensate. 

It is interesting to note that Beckertt & Stolper (1994) 
have found that by melting a CaO-MgO-Al 2 3 -Si0 2 -Ti0 2 sys- 
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urface reactions r assumed to form the solid materials s. The efficiency of the reaction is limited by the 
:ey species, which has the lowest abundance among the reactants. The notation % in the r.h.s. column means 
nd collision (and sticking) event initiates one reaction (see v^ ey in Eqs. |4] and ITOb . Data sources for the 
i (and saturation vapour pressures): (1) Helling & Woitke (2006); (2) Nuth & Ferguson (2006); (3) Sharp & 



Solid s 



Surface reaction 



Key species 



Ti0 2 [s] 
rutile 

(1) 



Ti0 2 — > Ti0 2 [s] 
Ti + 2 H 2 — > Ti0 2 [s] + 2H, 
TiO + H 2 — > Ti0 2 [s] + H 2 
TiS + 2 H 2 — > Ti0 2 [s] + H 2 S 



H 2 



Ti0 2 
Ti 
TiO 
TiS 



Si0 2 [s] 
silica 
(3) 



Si0 2 — > Si0 2 [s] 

SiO + H 2 — > SiO,[s] + H 2 

SiS + 2 H 2 — » Si0 2 [s] + H 2 S + H 2 



Si0 2 

SiO 

SiS 



SiO[s] 
silicon mono-oxide 
(2) 



SiO — -» SiO[s] 

Si0 2 + H 2 — » SiO[s] + H 2 

SiS + H 2 — > SiO[s] + H 2 S 



SiO 

Si0 2 

SiS 



Fe[s] 
solid iron 
(1) 



Fe — > Fe[s] 

FeO + H 2 — > Fe[s] + H 2 
FeS + H 2 — » Fe[s] + H 2 S 
Fe(OH) 2 + H 2 — » Fe[s] + 2 H 2 



Fe 

FeO 

FeS 

Fe(OH) 2 



FeO[s] 
iron (II) oxide 
(3) 



FeO — > FeO[s] 
Fe + H 2 — > FeO[s] + H 2 
FeS + H 2 — > FeO[s] + H 2 S 
Fe(OH) 2 — > FeO[s] + H 2 



FeO 

Fe 

FeS 

Fe(OH) 2 



FeS[s] 
iron sulphide 

(3) 



FeS — » FeS[s] 

Fe + H 2 S — > FeS[s] + H, 

FeO + H 2 S — > FeS[s] + H 2 

Fe(OH) 2 + H 2 S — > FeS[s] + 2 H 2 



FeS 
Fe 

minjFeO, H 2 S) 
min(Fe(OH) 2 , H 2 S) 



Fe 2 3 [s] 
iron (III) oxide 
(3) 



2 Fe + 3 H 2 —> Fe 2 3 [s] + 3 H 2 

2 FeO + H 2 — > Fe 2 3 [s] + H 2 

2 FeS + 3 H 2 — > Fe 2 3 [s] + 2 H 2 S + H 2 

2 Fe(OH) 2 — > Fe 2 3 [s] + H 2 + H 2 



>/ 2 Fe 
feFeO 
>/ 2 FeS 
>/ 2 Fe(OH) 2 



MgO[s] 
periclase 
(3) 



MgO — > MgO[s] 
Mg + H 2 — > MgO[s] + H 2 
2 MgOH — > 2 MgO[s] + H 2 
Mg(OH) 2 — > MgO[s] + H 2 



MgO 
Mg 

feMgOH 
Mg(OH) 2 



MgSi0 3 [s] 
enstatite 

(3) 



Mg + SiO + 2 H 2 —> MgSi0 3 [s] + H 2 

Mg + SiS + 3 H 2 — > MgSi0 3 [s] + H 2 S + 2 H 2 

2 MgOH + 2 SiO + 2 H 2 — > 2 MgSi0 3 [s] + 3 H 2 

2 MgOH + 2 SiS + 2 H 2 — > 2 MgSi0 3 [s] + 2 H 2 S + 2 H 2 

Mg(OH) 2 + SiO — > 2 MgSi0 3 [s] + H 2 

Mg(OH) 2 + SiS + H 2 —> MgSiQ 3 [s] + H 2 S+ H 2 



minjMg, SiO) 
minjMg, SiS) 
minjfeMgOH, feSiO) 
minjfeMgOH, feSiS) 
min(Mg(OH) 2 , SiO) 
min{Mg(OH) 2 , SiS) 



Mg 2 Si0 4 [s] 
forsterite 

(3) 



2 Mg + SiO + 3 H 2 — > Mg 2 Si0 4 [s] + 3 H 2 
2 MgOH + SiO + H 2 — > Mg 2 Si0 4 [s] + 2 H 2 
2 Mg(OH) 2 + SiO — > Mg 2 Si0 4 [s] + H 2 + H 2 
2 Mg + SiS + 4 H 2 — > Mg 2 Si0 4 [s] + H 2 S + 3 H 2 
2 MgOH + SiS + 2 H 2 — > Mg 2 Si0 4 [s] + H 2 S + 2 H 2 
2 Mg(OH) 2 + SiS — > Mg 2 SiQ 4 [s] + H 2 + H 2 S 



minjfeMg, SiO) 
minj'/ 2 MgOH, SiO) 
minj'/ 2 Mg(OH) 2 , SiO) 
minjfeMg, SiS) 
min{'/ 2 MgOH, SiS) 
min{V 2 Mg(OH) 2 , SiS) 



Al 2 3 [s] 
aluminia 

(3) 



2 Al + 3 H 2 — 
2 AlOH + H 2 • 
2 A1H + 3 H 2 • 
A1 2 + 2 H 2 - 
2 A1S + 3 H 2 - 



Al 2 3 [s] + 3H 2 
-> Al 2 3 [s] + 2H 2 
-» Al 2 3 [s] +4H 2 
-> Al 2 3 [s] + 2H 2 
-» Al 2 3 [s] + 2H 2 S + H 2 



2 A10 2 H Al 2 3 [s] + H 2 



y 2 Al 
fcAlOH 

V 2 aih 

A1 2 
feAlS 
'/ 2 A10 2 H 



CaTi0 3 [s] 
perovskite 

(3) 



Ca + TiO + 2 H,0 
Ca + Ti0 2 + H 2 - 
Ca + Ti + 3 H 2 — 
CaO + Ti + 2 H 2 • 
CaO + TiO + H.O - 



^ CaTi0 3 [s] + 2H 2 
^CaTi0 3 [s] +H 2 
CaTi0 3 [s] + 3H 2 
->CaTi0 3 [s] + 2H 2 
^CaTi0 3 [s] +H 2 



CaO + Ti0 2 —> CaTi0 3 [s] 
CaS + Ti + 3 H 2 CaTi0 3 [s] + H 2 S + H 2 
CaS + TiO + 2 H 2 CaTi0 3 [s] + H,S + 2 H 2 
CaS + Ti0 2 + H 2 — > CaTi0 3 [s] + H 2 S 



minjCa, TiO) 
minjCa, Ti0 2 ) 
minjCa, Ti) 
minjCaO, Ti) 
minjCaO, TiO) 
minjCaO, Ti0 2 ) 
minjCaS, Ti) 
minjCaS, TiO) 
minjCaS, Ti0 2 ) 
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tern, complex compounds like CaA^O?, Ca3Ti2Ai2Si20i4 and 
Ca3Al2Si40i4 can form. It seems logic to conclude that sim- 
pler compounds like MgO etc. need to form before more com- 
plex compounds can be built. It might furthermore be interest- 
ing to consider Earth-crust-like solids like CaSiOi or CaCC>3. 
CaSiCb is very abundant on Earth but it forms only under high- 
pressure conditions inside the Earth crust at 10 4 - 10 6 bar (A. P. 
Rossi 2006, priv. com.). CaCCb (carbonaceous calcite), as all 
the other carbon bearing solids solids, can only form in an 
oxygen-rich environment if the carbon is not locked into the 
CO molecule. Under the high pressure conditions in brown 
dwarf atmospheres, this release takes place only below w 800K 
under chemical equilibrium conditions. Toppani et al. (2005) 
have concluded that carbonates can only form in a H20(g)- 
CC<2(g)-rich, high-density region under conditions far away 
from thermal equilibrium. 

Our selection of 12 solids is assumed to be formed by 60 
chemical surface reactions (see Table [T]). We assume for the 
sticking coefficient a,- = 1 due to the lack of data (see discus- 
sion in Paper V). The Gibbs free energy data of the solid com- 
pounds, from which the supersaturation ratios S are calculated, 
is mostly taken from Sharp & Huebner (1990) which were ob- 
tained for crystalline materials. The data for Ti02[s] is given in 
Paper III. The new vapour pressure data for SiO (Ferguson & 
Nuth 2006) is used. The reaction supersaturation ratios S r are 
calculated as outlined in Appendix B of Paper V. 

4. Results 

4.1. Cloud formation and cloud structure 

The vertical structure of the dust cloud layer in a brown dwarf 
atmosphere is depicted in Fig. Q] The results are shown for a 
stellar parameter combination typically associated with the L 
dwarf regime. 

The cloud structure results from a hierarchical dominance 
of nucleation (uppermost layers), growth & drift (intermediate 
layers), and evaporation (deepest layers). In the uppermost lay- 
ers, small grains of size ~ 0.01 fim form by nucleation and grow 
further as they settle down the atmosphere. The particles reach 
a maximum size of about 200 jjm at cloud base, shortly before 
they evaporate completely. The grain fall speeds first decrease 
with increasing depth because of the increasing ambient gas 
densities, and then re-increase as the particles grow rapidly. 
Finally, they fall faster than they can grow (rain) and reach 
an approximately constant fall speed of a few m/s. Eventually, 
the grains enter the hotter atmospheric layers where they are 
no longer thermally stable. Hence, the grains shrink in size 
and dissolve into the surrounding hot and convective gas (the 
Schwarzschild pressure, where the atmosphere becomes con- 
vectively unstable, is 3.54bar in this model). These results re- 
flect the stationary character of the dust component in sub- 
stellar atmospheres, where dusty material constantly forms at 
high altitudes and settles downward, and simultaneously, fresh 
uncondensed material is mixed up by convective motions and 
overshoot (see Ludwig et al. 2002, 2006; Young et al. 2003). 
These general results resemble well the results of Paper III, 
where only one sample dust species Ti02[s] was considered. 
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Fig. 1. Calculated cloud structure for a model with T e g = 
1800K and log g = 5 (cond AMES). 1 st panel: prescribed 
temperature T (solid) and mixing time scale r m ; x (dashed). 2 nd 
panel: nucleation rate 7+ (solid) and net growth velocity Xnet 
(dashed). 3 rd panel: mean grain size (a) - -\J3 / '(4-n) L\ j 'Lq 
(solid) and mean fall speed (v dr ) = yjngpdia) / (2pcj) 
(dashed). The 4 th panel depicts the effective supersaturation ra- 
tios S e ff for all solids involved. 



However, in comparison to Paper III, there are new fea- 
tures resulting from the inclusion of more than one solid growth 
species and the consideration of the more abundant conden- 
sates. The net growth speed of the grains Xmt 1S never deter- 
mined by a single species alone, but results from a complicated 
superposition where the individual contributions ^ et can be 
positive (growth) or negative (evaporation). The little kinks in 
Xnet (see 2 nd panel in Fig. [TJ result from the evaporation of one 
material which becomes thermally unstable at a certain tem- 
perature, in this case Mg2Si04[s] at around 1800 K and Fe[s] 
at around 2000 K. 

4.2. Vapour saturation in the atmosphere 

In order to discuss deviations from phase-equilibrium, we con- 
sider the supersaturation ratio S , which usually indicates net 
growth for S > 1 and net evaporation for S < 1 . However, for the 
complicated growth reactions listed in Table [T] where mostly 
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Fig. 2. Volume fractions V s /V tot (top) and mass fractions 
M s /M tot (bottom) of the various solid compounds for the model 
depicted in Fig. Q] The temperature scale follows the atmo- 
spheric stratification being cool at the top and hot at the bottom. 

at least 2 gas particles need to collide with the grain's surface 
(called "type III reactions" in Paper V), the supersaturation ra- 
tio is not unique, but reaction-dependent (see Appendix B of 
Paper V). In addition, the Z? s s urf -factors are involved to account 
for the inequality of the active surfaces for growth and evapo- 
ration (see Paper V). The thermal stability of the solids is then 
better defined by^ et = (see Eq. |4j». In order to discuss the sat- 
uration of the atmosphere anyway, we define a unique effective 
supersaturation ratio S e ff for each solid as 
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These effective supersaturation ratios S e ff shown in Figure Q] 
(4 th panel) demonstrate that actually none of the considered 
solids obeys the condition of phase-equilibrium (S e ff = 1) 
throughout the entire cloud layer in our model. In particular, 
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Fig. 3. Volume fractions V s /Vtot for different T e g (top) and 
different log g (bottom). Only the more abundant solid com- 
pounds are viewed in contrast to Fig. [2] 

the atmospheric gas is strongly supersaturated in the uppermost 
layers where the nucleation takes place. However, in contrast 
to TiC>2[s] discussed in Paper III, the more abundant silicates 
and iron compounds reach a state of quasi-phase-equilibrium 
(S e ff ^ 1) already at p = 0.05 ... lbar whereas the rare conden- 
sates like Ti02[s], A^Chts], and CaTiC>3[s] practically never 
obey the phase-equilibrium condition. This behaviour is caused 
by the different element depletion time scales being longest for 
the low abundant elements as we have demonstrated in Paper V 
(Sec. 4.4.3.). Since the observable molecular features are typ- 
ically formed high in the atmosphere and even the broad dust 
features in the wavelength region 7-20 pm originate from lay- 
ers p x 0.2 ... 0.6 bar in this model (see Sect. 15.2b we conclude 
that the observable dust in brown dwarf atmospheres in not in 
phase-equilibrium with the gas. 

4.3. Chemical composition of the cloud particles 

The chemical composition of the dust grains is expressed by the 
solid material volume fractions Vs/Vtot (Fig- El top panel) and 



8 



Ch. Helling et al.: Dust in Brown Dwarfs and Extra-solar Planets 



Table 2. Typical dust volume composition V s /V tot , mass com- 

position M s /M tot , and mean particles sizes (a) [fim] as function °g 

of local temperature for models as depicted in Fig. [2] (/*) - \ 

increasing in T-interval; (\) - decreasing in T-interval. ^ 
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MJV tBt 
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24% MgSi0 3 , 20% Mg 2 Si0 4 


l.MgSi0 3 


10 i 
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(★) Mean grain sizes {a} can reach up to 10 3 5 in lower gravity atmo- 
spheres like giant-gas planets (compare Fig.f4]l. 
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mass fractions M s /M tot = V s Ps/(VtotPd) (l° wer panel). The dust 
material composition is mainly controlled by the temperature, 
and shows a constantly re-occurring pattern in models with dif- 
ferent stellar parameter for brown dwarfs and gas giant-planets 
(see also Table[2]i. Once a solid becomes thermally unstable, the 
respective elements evaporate into the gas, making them avail- 
able again to form other, more stable condensates. In this way, 
we find a complicated mix of all condensates high in the atmo- 
sphere with volume fractions resulting from kinetic constrains 
that favour the formation of the more abundant silicates and Fe- 
oxides. As these dirty particles settle down the atmosphere they 
stepwise purify, until only the most stable component parts like 
Ai2(D3[s], Fe[s], and CaTi03[s] remain. 

Therefore, two classes of solids can be distinguished: the 
high-temperature condensates Fe[s] and AhCbts] with some 
contributions of Ca-Ti-oxides in the deeper layers, and the 
medium-temperature condensates MgSK^s], Mg 2 SiC>4[s] 
and SiC>2[s] with some SiO[s] and FeS[s] in the upper layers. 

As the temperature increases, SiO[s] starts to evaporate and 
SiChts] becomes the second most abundant solid material by 
volume fraction. Next, MgSiCbts] evaporates which sets free 
some Mg and Si to form more Mg2Si04[s] and SiC>2[s], mak- 
ing SiC>2[s] again the second most abundant solid. The tran- 
sition from the medium-temperature to the high-temperature 
composition is characterised by the evaporation of Mg2Si04[s], 
which leaves only Fe[s], A^Chts], TiC>2[s] and CaTi03[s] as 
stable condensates. As Fe[s] evaporates, the cloud base is soon 
reached where eventually the remaining Al-Ca-Ti oxides evap- 
orate. 



Fig. 4. Comparison of results for different stellar parameter. 1 st 
box: mass density p, 2 nd box: nucleation rate /+, 3 rd box: net 
growth velocity x aet 4 th box: mean particle size (a), 5 th box: 
mean fall speed (v dr ). All quantities are plotted as function of 
temperature T. Dotted line: T e g - 1300 K, log g = 5; solid line: 
r eff = 1800K,log g = 5; dash-dotted line: r eff = 1300K, log g = 
3 ; dashed line: T eS = 1 800 K, log g-3. 

The volume fractions are important for the dust opacities 
(see Sect. |5), and in the observable layers the Mg-silicates turn 
out to be the most relevant condensates. For completeness we 
note, however, that the mass fraction M s /M tot is dominated by 
Fe[s] from about 1200 K to 2000 K. 

We recognise that the chemical composition of the dust 
grains does to some extend depend on the completeness of the 
selection of solids and on the availability of material quanti- 
ties (like e.g. the sticking coefficients as shown in Sect. 4.6. 
in Paper V). For example, the partial volume of Mg2SiC<4[s] 
and MgSiC<3[s] decrease if MgO[s] is included, and the partial 
volume of SiC<2[s] decreases if SiO[s] is included. For this rea- 
son, the partial volume of SiC<2[s] is smaller than described in 
Helling et al. (2006). 

4.4. Dependence on T eS and log g 

Figure [3] shows the dependence of the dust material composi- 
tion on r e ff and log g. If plotted against temperature, the vol- 
ume fractions show a robust pattern for all calculated models 
(see Table|2|i. Only a slight shift to higher temperatures (deeper 
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layers) can be noticed for lower r e g (upper plot) and higher 
log g (lower plot). 

An increase of log g results in a more compact atmosphere, 
i. e. generally higher pressures in the atmosphere. Since the 
dust sublimation temperatures increase with increasing pres- 
sure, the dust remains stable to even higher temperatures for 
larger log g. For the same reason, models with lower T e $ 
show dust at comparably higher temperatures. For lower T e s, 
a certain temperature is reached deeper inside the atmosphere, 
where the pressure is higher and, hence, the dust is more stable. 

Figure |4] shows the dependence of the nucleation rate, the 
mean particle size and the mean fall speed on T e g and log g. 
The nucleation zone lies typically between 600 K and 1400 K 
in all models. The maximum rate reaches higher values for 
higher log g because of the higher gas densities. For lower 
log g, the nucleation maximum is more extended and shows 
a more complicated shape that is probably caused by the closer 
neighbourhood to the convective zone, which makes the up- 
mixing of fresh elements for nucleation more likely. 

Concerning the mean particle size, all models show about 
the same small particles of order 0.01 fim high in the atmo- 
sphere, which grow to large particles between about 100 /j and 
1000 fi in the deep layers. The maximum particle sizes reached 
at cloud base are larger for lower log g. As the nucleation zone 
ends, there is a zone of rapid grain growth around 1400 K to 
1700 K. In this zone, most of the solid particles are actually 
"mixed away" according to our simple approach to treat the 
mixing by convective motions and overshoot. The small num- 
ber of particles that stay, however, settle down deeper into a 
denser environment where elements are mixed up with higher 
efficiency due to the decreasing distance from the convection 
zone. The models show that this small number of growing par- 
ticles is sufficient to maintain a state close to phase equilibrium 
with the gas concerning most elements, i. e. the growth is ex- 
haustive. Thereby, the particles grow further along their way 
down the atmosphere while their number is ever decreasing due 
to lethal mixing. 

The growth of the particles, however, is limited by their 
own fall speed which increases with size. The grains eventu- 
ally reach a size where their residence time is so small that 
further growth becomes negligible. The residence time scale 
Tsink = Vfc/Hp, however, depends on the atmospheric scale 
height Hp which is 100 times larger for the log g = 3 mod- 
els. Therefore, the dust particles have more time to grow in the 
low log g models, producing larger sizes and larger fall speeds. 

We note that the Knudsen numbers fall short of unity in 
the deeper layers, for particles larger than about 1 fim in the 
log g-5 models and about 10 fim in the log g = 3 models. This 
means, that the frictional force and the growth velocity should 
be calculated for the case of small Knudsen numbers in the 
deeper layers, making necessary a Knudsen number fall differ- 
entiation (see Paper II for details). Therefore, the results of this 
paper concerning the particle sizes in the deeper layers must be 
taken with care. The true fall speeds are probably larger there 
and the true growth velocities are probably smaller as com- 
pared to our results, i. e. we expect that (a) remains smaller in 
the deeper layers as compared to Fig.|4] 
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Fig. 5. Dust-to-gas mass ratios pd/Pgas and dust mass column 
densities Md lls t for brown dwarf and giant-planets atmospheres. 
1 st panel: linear plot, 2 nd panel: logarithmic plot. Dotted line: 
r eff = 1300K, log g = 5; solid line: T eS = 1800K, log g = 5; 
dash-dotted line: T e g = 1300 K, log g = 3; dashed line: T e g = 
1800 K, log g=3. 

Table 3. Maximum dust-to-gas ratio and temperature interval 
where the dust-to-gas ratio is larger than half maximum. 



T* [K] 


log g 


(— ) [io- 3 ] 

^Puus 'max 


temp, interval [K] 


1300 
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2.6 


1300 - 1800 


1800 


5 


2.6 


1300 - 1800 


2200 


5 


2.1 


1450 - 1750 


2500 


5 


0.2 


1550 - 1800 


1300 
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4.0 


950 - 1600 


1800 


3 


3.8 


1000 - 1600 


2200 
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2.3 


1300 - 1500 


2500 


3 


0.04 


1550 - 1700 



4.5. Dust-to-gas ratio 

Figure|5]shows the dust-to-gas mass ratios pd/Pgas and the dust 
mass column densities M<just = / P&dz [g/cm 2 ] as function of 
pressure. The dust-to-gas ratio generally first increases inward 
log-linear, then reaches a plateau and finally decreases rapidly 
as the dust becomes thermally unstable. The upper plot shows 
the same quantities on a linear scale, more emphasising the 
truly dusty layers. 

For the four models depicted in the previous figures, the 
dust-to-gas ratio reaches a constant maximum value between 
about 0.25% and 0.4% in a temperature interval T e [(1000 - 
1300) K... (1600 - 1800) K] rather independent of T eff , where 
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the temperature interval boundaries increase with increasing 
log g (see Sect. |4.4| i. The width of the dusty temperature win- 
dow is about 500 K to 600 K for the depicted models with 
r eff = 1300 K and 1800 K (compare Table EJ. We note that this 
result is actually quite consistent with the simple dust approach 
made by Tsuji (2002), although the window is broader than as- 
sumed by Tsuji according to our findings. Similar maximum 
values are found even if we increase T e g to 2200 K. However, 
if we increase T e s further to 2500 K, the maximum dust-to-gas 
ratio drops by more than one order of magnitude (see Tabled. 
We conclude that T e ff » 2200 K is the threshold value for truly 
dust-rich layers to occur in the models presented here. 

The resulting dust-to-gas ratios demonstrate that the cloud 
layer is primarily attached to the local temperature. For lower 
J e (f, the dust layer sinks in deeper into the atmosphere and 
eventually disappears from the observable layers. Since the 
dust is then present in denser regions, the dust mass column 
density Md us t increases. 

4.6. Remaining gas-phase chemistry & metallicity 

The remaining gas-phase composition depends on the amount 
of elements not locked up into dust grains. As a result of our 
model, phase-equilibrium is not valid in the upper layers (see 
Fig-El 4 th panel) and should not be used to determine gas parti- 
cle abundances. The remaining gas abundances e; are strongly 
sub-solar in an extended layer above the cloud layer where the 
pressure drops by about 3 orders of magnitude (about 10 scale 
heights), see Figs. [6] and [7] Inside the cloud layer, the metal 
abundances increase with increasing atmospheric depth and fi- 
nally reach even slightly larger than solar values at cloud base, 
where those elements are released by evaporation that have 
been locked up into grains in the upper layers. 

Figure [6] shows this phase lag between metal gas abun- 
dances and dust abundances clearly. Considering a path from 
the top to the bottom of the atmosphere, first the metals disap- 
pear, then the dust appears, then the metals reappear and finally 
the dust disappears (the sum of dust and gas abundances is not 
constant). 

The calculation of the gas phase element abundances al- 
low for a discussion of the metallicity in dust forming atmo- 
spheres (see Fig. O, i. e. the logarithmic differences of the gas 
abundances with respect to the solar values. The metallicities 
are not the same - not even similar - for different elements. 
The strongest depletions occur for Ti, Mg and Si, for which 
the metallicities reach values as low as -6 to -7. The deple- 
tion of S and Ca is much less significant (> -0.5) in com- 
parison. Concerning Ca, this is questionable because it is due 
to our limited choice of solids. There is only one Ca bearing 
solid species considered, CaTi3[s], and since Ti is less abun- 
dant than Ca, Ca cannot be locked up completely into grains. It 
is important to note, however, that according to our models, the 
gas depletion factors are kinetically limited by a maximum fac- 
tor of about 10 6 . Larger depletion factors simply cannot occur 
because the depletion timescale (see Paper V) would exceed 
the mixing timescale. In contrast, complete phase equilibrium 
would imply that the depletion factors can reach values as high 
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Fig. 6. Element abundances in dust ea (thick solid) and in the 
gas phase e\ (dashed). The solar values e^sun (thin solid) and 
the linear dust-to-gas ratio pd/Pgas (lowest panel) are shown for 
comparison. The stellar parameter are T e g = 1 800 K and log g = 
5 like in Fig. Q] 



as 10 50 for low temperatures. Figure [7] also demonstrates that 
the metallicities change with altitude and depend on the stellar 
parameter. 

A special feature of our models is that the metal abundances 
in the gas phase re-increase high above the cloud layer. In our 
model, all elements are constantly mixed upward, and although 
the mixing timescale is as long as T m [ x = 10 10 s in the high- 
est layers, the timescale for a gaseous particle to find a surface 
to condense on exceeds this timescale, because the dust abun- 
dance drops steeply above the cloud layer. Consequently, the 
metal abundances asymptotically re-approach the solar values 
high above the cloud layer. The only exception is Ti, which 
can disappear from the gas phase via nucleation. This fea- 
ture distinguishes our models from all other published models, 
and possibly opens up a possibility to discriminate between 
the models by detailed line profile observations of neutral al- 
kali and alkaline earth metal atoms. Our dust treatment results 
in deeper and broader resonance lines of Na I, and K I in the 
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not presentable anymore. 



red part of the spectrum, which form high in the atmosphere 
(Johnas et al. 2007). 

4. 7. Grain size distribution 

In Fig. [8] we show the results for the potential exponential grain 
size distribution function [cirT 4 ] (see Appendix IA.2I ) 

with parameters derived from the calculated dust moments 
Lj(z). A number of altitudes has been selected arbitrarily 
throughout the entire dust cloud for visualisation. 

The evolution of the grain size distribution through the 
cloud layer are similar for models with the same log g. 
Considering the log g — 5 models (upper row in Fig [8j, the 
size distribution starts relatively broad in the upper atmosphere 
where the nucleation is active, because the constant creation 
and simultaneous growth of the particles causes a broad dis- 
tribution. Once the nucleation ceases, however, all particles 
merely shift in size space by a constant offset Aa due to fur- 
ther growth, which means a narrowing in A(loga). Since the 
particles grow by more than 4 orders of magnitude on their 
way down the atmosphere, the size distribution finally becomes 
strongly peaked. 



In the log g = 3 models (lower low in Fig [8), the evolution 
of the size distribution function is more complicated. As shown 
in Fig.|H the nucleation rate has a a small shoulder as function 
of pressure on the left hand side. This feature leads to a narrow- 
ing of f(a) at first, followed by a re-widening of the size dis- 
tribution function with increasing depth. The nucleation zone 
is closer to the convective zone and sometimes even overlaps 
with the convective zone in the log g — 3 models. Therefore, al- 
though the nucleation rate becomes tiny with increasing depth, 
it does not vanish completely until about 1600 K which still 
influences the size distribution and keeps it broad. 



5. Simple radiative transfer modelling 

5.1. Model description 

We simulate the grain absorption features in two brown dwarf 
atmosphere cases (T e g = 1800K, log g — 5 and = 1300K, 
log g = 5) and one example for a giant-gas planet atmosphere 
(Teff - 1300K, log g — 3) using a simple radiative transfer 
code. The code considers the solid-state opacities similar to 
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Fig. 7. Element abundances involved in the dust formation pro- 
cess for two effective temperatures and two different gravities, 
dotted blue: r eff = 1300 K, log g = 5; solid green: r eff = 1800K, 
log g — 5; dash-dotted blue: - 1300 K, log g — 3; dashed 
green: T eS = 1 800 K, log g = 3 . 

Helling et al. (2006fl Grain opacities are calculated accord- 
ing to their size distribution f(a,z) [cm -4 ] and their solid ma- 
terial volume composition V s (z) (Figs. [2] [8} using the effec- 
tive medium theory according to Bruggeman (1935) to calcu- 
late the effective optical constants and Mie theory for spheri- 
cal particles to calculate the extinction efficiencies Q ex t(a, A). 
The main difference between the previous code and the present 
one is that the grains are here assumed to be distributed ac- 
cording to the potential exponential size distribution function 
(see Appendix IA.21 > determined from the calculated dust mo- 
ments at height z, whereas in the previous model a (^-function, 
f(a,z) = pL (z)5(a - (a)(z)) representing the mean particle 
size (a), was assumed. The present code includes the 12 solid 
species discussed in the previous sections (see Table l5.2l for the 



3 More complex methods applicable to hydrodynamic environ- 
ments are given in e.g. Helling & J0rgensen (1998) and Woitke (2006). 



references for the optical constants). The total grain extinction 
cross-section (see inner integral in Eq. [T3l is computed using 
a Gauss-Legendre quadrature integration over the entire grain 
size distribution function at given atmospheric height. 

We focus the radiation transfer on the 7/vm to 20^/m wave- 
length range, which encompasses the Si-O stretching mode 
(centred at 9.7 fim) and Si-O bending mode (around 18 //m) of 
crystalline silicates (Mg2Si04 and MgSiOs) and quartz (SiOi). 
The model calculates the optical depth of the dust component 



dust 



(Z) 



Z °o 

// 



f(a,z') na 2 Q ext (a,A,V s (z')) dadz' (13) 



and determines the geometrical depth zo(^) where rj| ust (zo) = 1 
for each wavelength. Furthermore, we replace the opacity of 
the solid species by vacuum to estimate their respective effects 
on the output spectrum. To calculate the transmission spectra 
(first row Fig. |5), we remove a blackbody emission as contin- 
uum (Z?i(rbb)) from the spectrum (F A ) as 

transmission = F A (r[zoW)]) /fi^bb) , (14) 

i. e. all transmission curves for one model atmosphere are di- 
vided by the same black body continuum: (T e ff, log g; Tbb) = 
{(1800K, 5.0; 1310K), (1300K, 5.0; 1205K), (1300K, 3.0; 
910K)}. 



5.2. Results 

The absorption features can be distinguished from the contin- 
uum in all cases. But they amount to a maximum of 6% only 
at 9.7yum in all three cases. The positions of the absorption fea- 
tures (9.7yizm and 17-18 /mi) are typical of pure absorption by 
amorphous silicates, which implies that scattering by the larger 
grains in the size distribution does not contribute significantly 
to the total cross-sections. This result differs from the single 
grain size situation where scattering may modify the features, 
e. g. by shifting the peak absorptions. 

The temperature and pressure at zq(A), where rj| ust = 1, 
are plotted below the transmission spectra. The weak absorp- 
tion reflects the shallow temperature and pressure variations in 
the r dust =s 1 region. The contribution of quartz is negligible, 
as testified by the same output spectrum whether quartz opac- 
ity (dash-dot line) is present or not. The two silicate species 
MgSi03 and Mg2Si04 contribute about equally to the output 
spectrum (dotted and dashed line). 

Although grains exist even high in the atmosphere in this 
model, their number densities are too small to effect the opac- 
ity until the region of rapid grain growth is reached. The 
rapid grain growth implies a concomitant sudden rise in opac- 
ities. The phenomenon is particularly acute in the (T e ff = 
1800K, log g = 5) model. Therefore, the absorption features 
probe in particular the cloud deck. The inclusion of gas opac- 
ities in the radiative transfer model will complicate the 7-20 
fim spectra, definitively rendering the search for silicates fea- 
tures in brown dwarf atmospheres difficult. 
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Table 4. Reference for the optical constant of amorphous ma- 
terials used in the radiative transfer modelling. 



Solid species 


Reference 


Ti0 2 [s] 


Ribarsky et al. (1985) 


Q\C\ Tel 


rienning et al. yvyyi ) 


SiO[s] 


Philipp in Palik(1985) 


Mg 2 Si0 4 [s] 


Jager et al. (2003) 


MgSi0 3 [s] 


Dorschner et al. (1995) 


MgO[s] 


Hofmeister et al. (2003) 


Al 2 3 [s] 


Begemann et al. (1997) 


Fe[s] 


Ordaletal. (1985) 


FeO[s] 


Henning et al. (1995) 


FeS[s] 


Begemann et al (1994) 


Fe 2 3 


http ://w ww. astro . uni-j ena.de/Laboratory/ 




OCDB/index.html 


CaTi0 3 [s] 


Posch et al. (2003) 



6. Conclusive summary 

Based on a detailed micro-physical description of the for- 
mation of composite dust particles including nucleation, 
growth/evaporation and drift we have investigated the forma- 
tion, chemical composition and spectral appearance of quasi- 
static cloud layers in quasi-static brown dwarf and gas-giant 
atmospheres. Our models show that 

- There is one single cloud layer which stretches from about 
(1000- 1300) K to (1600-1800) K, which shifts to slightly 
higher temperatures for higher log g values. 

- The material composition of the cloud particles changes 
continuously through the cloud layer. The seed particles 
created high in the atmosphere are simultaneously covered 
by all kinds of solid materials, favouring the abundant con- 
densates like MgSi0 3 [s], Mg 2 Si0 4 [s], Si0 2 [s] and Fe[s]. 
As the particles sink in deeper, the higher temperatures en- 
force partial evaporation which stepwise purifies the grains, 
leaving only the most stable condensates like Fe[s] and 
Ai2C>3[s] at cloud base. 

- The material composition described above is primarily a 
function of local temperature and robust against changes in 
r eff and log g. 

- The mean sizes of the cloud particles increase continuously 
through the cloud layer. The particles are found to be very 
small in the high atmospheric layers ((a) * 0.01 /mi) with 
a relatively broad size distribution, but as large as several 
100 /mi with a strongly peaked size distribution at cloud 
base, where they finally shrink due to evaporation and dis- 
solve into the surrounding gas. The maximum particle sizes 
reached are larger for smaller log g. 

- The gas is highly supersaturated in the upper atmosphere, 
where the nucleation takes place, but reaches a state close 
to saturation from the cloud deck downward concerning the 
abundant medium-temperature condensates. For high alti- 
tudes in general and for high-temperature condensates in 
particular, phase equilibrium is not valid, since the deple- 
tion timescale exceeds the mixing timescale. 



- No cloud layer will form without the formation of seed par- 
ticles, which is different from terrestrial planets where seed 
particles can be swapped up from the crust. 

- The consequence of the nucleation process is that the seed 
forming elements are stronger depleted from the gas phase 
in the high atmosphere than those elements which need 
an alien surface to condense on, which leads to a highly 
Ti-depleted upper atmosphere in our model. This finding 
could be useful to identify first nucleation species in brown 
dwarfs and gas-giants. 

- Our models predict a much smaller degree of gas phase 
depletion as compared to phase-equilibrium models. The 
maximum metallicity depression is about 6 orders of mag- 
nitude in our model for Mg, Si, Fe and Al. Therefore, age 
estimates based on phase-equilibrium models may have 
overestimated the brown dwarf's age considerably. 

- As a result of the phase-non-equilibrium, the carriers of the 
metallic resonance lines, e.g. Nal, KI, remain more abun- 
dant at high altitudes which should make these lines deeper 
and broader in the red part of the spectrum. This has by now 
been shown in Johnas et al. (2008). 

- The oxygen-depletion is not large enough to turn the pre- 
mordially oxygen-rich substellar atmosphere into a carbon- 
rich atmosphere with a carbon-to-oxygen ratio > 1 . 

- The broad dust absorption features attributed to the Si-O 
stretching mode (centred at 9.7 /mi) and the Si-O bend- 
ing mode (around 18//m) should be present in brown 
dwarf spectra, but their maximum absorption is weak (< 
6%) and hence difficult to detect. The spectra primarily 
probe the cloud deck between about 850 K and 1300K, 
where Mg2Si04[s] and MgSi03[s] are the most important 
solid species. The positions of the dust absorption features 
(9.7/mi and 17-18/mi) are typical of pure absorption by 
amorphous silicates, which implies that scattering by the 
larger grains in the size distribution does not contribute sig- 
nificantly to the total cross-sections. 



Appendix A: Grain size distribution function 

In order to find the dust particle size distribution function 
f(a) [ctrT 4 ] we switch to another set of dust moments Kj 
where the integrals are performed in radius-space a (Dominik 
et al. 1986, Gauger et al. 1990) rather than in volume space V 
(Dominik et al. 1993). Substitution for V-(4na 3 )/3 yields 

Kj = J f(a)a J da = ^—j pLj (A.l) 

As argued in the main text, we are left with four known proper- 
ties of the size distribution function, namely K l , K2, K3 and K4. 
We will introduce a suitable functional formula for f(a) with a 
set of four free coefficients and will determine these coefficients 
from the known dust moments. This exercise is carried out for 
two functions in the following. 
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A .1. Double delta-peaked size distribution function 

One option is to consider the superposition of two Dirac- 
functions 

f(a) = Nt 6(a - a\) + N2 S(a - CI2) , (A.2) 

where 5 is the Dirac-function, N\,N2 [cm -3 ] are two dust par- 
ticle densities and a\,a% [cm] are two particle radii. 

Assuming a\ > at and a 2 > at, where at is the lower in- 
tegration size corresponding to Vt, the moments are given by 
Kj = Nia^ + N2aL from which we can determine the four free 
coefficients. a\ is found to be the positive root of 

a 2 1 (Kl-KiK 3 ) + ai(K 1 K 4 -K2K 3 ) + (Kl-K 2 K 4 ) = (A.3) 

and the other free coefficients are given by 
a\K 2 - K 3 



a 2 = 



Ni = 



CL\Kl — K2 



(aiKi - K 2 y 



{a x K 2 -K 3 ){K 3 -2a x K 2 +a\K x ) 



N 2 



ai(K 3 -2aiK 2 +alKi) 



(A.4) 
(A.5) 
(A.6) 



A.2. Potential exponential size distribution function 

Another, more continuous option is to consider the function 

/(a) = a B exp(A-Ca) , (A.7) 

which, for positive coefficients A, B and C, is strictly positive 
with a maximum at B/C. Since there are only three coefficients 
in Eq. ( IA.7t we will determine them from K\, K 2 and ^3. 

To find a simple analytical solution, we extend the integra- 
tion in Eq. ( IA.1I ) to the interval [-00 ... +00] which usually in- 
troduces only a small error, since f(d) is quickly vanishing for 
small a (see Fig. 8). The result is Kj = exp(A + In F(B + 1 +j) - 
(B+1+ j) InC), by which the free coefficients can be deduced: 



B = 



2Ki K 3 - 3Kl 



C = (2+B)-^- 
K 2 

A = ln^i +(2+B)lnC-lnr(2+fi) 



(A.8) 

(A.9) 
(A. 10) 



where F is the generalised factorial function with F(n+ 1) = n\ 
for n € N and Y(x+ 1) = x V(x) for x e R + . 
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